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We present an evaluation of box diagram for the elastic ep scattering with proton in the intermedi- 
ate state. Using analytic properties of the proton form factors we express the amplitude via twofold 
integral, which involves the form factors in the space-like region only. Therefore experimentally 
measured form factors can be used in the calculations directly. The numerical calculation is done 
with the form factors extracted by Rosenbluth, as well as by polarization transfer methods. The 
dependence of the results on the form factor choice is small for Q 2 < 6 GeV 2 , but becomes sizable 
at higher Q 2 . 
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I. INTRODUCTION 

Study of the clastic electron-proton scattering 



is an important source of information about the proton structure. In the first order of the perturbation theory (PT) 
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' or Born approximation the reaction amplitude 

)£ ; M B om = —u^u ■ UT»(q)U (2) 

■ is expressed in terms of Dirac and Pauli proton form factors (FFs), Fi and F2, 

F^q) = F^q 2 )^ - F 2 (q 2 )[^, 7v ]£j, (3) 

^ I where u, U, u , U are 4-spinors of incoming and outgoing particles, a « 1/137 is fine structure constant, q is the 
momentum transfer from the electron to the proton. The linear combinations, Ge = F\ + 4^73-^2 and Gm = -Pi + F2, 
called the electric and magnetic FFs, are also widely used. 

During many years it was a common practice to extract FFs by the Rosenbluth separation method. FFs obtained by 
this method obey with good accuracy Ge/Gm = const for < Q 2 < 6 GeV 2 (Q 2 = —q 2 )- However since 2000 series 
of experiments was done using an alternative, polarization transfer method [J. These experiments yielded strikingly 
different results: Ge/Gm ratio decreased linearly with Q 2 . Since both methods are based on the Born approximation, 
the reason for discrepancy is likely the neglected higher order PT terms. There are two types of second order Feynman 
diagrams: first, diagrams, involving an exchange of only one virtual photon (vacuum polarization or vertex corrections) 
and second, two photon exchange or box diagram, Fig^a). The part of amplitude, coming from one-photon exchange 
diagrams, has the structure similar to 10) (certainly with another functions in place of Fi^) and cannot lead to the 
discrepancy between two methods. Such diagrams are usually taken into account in experimental analyses. 

Therefore the two photon exchange diagram plays the key role in understanding the experimental results and 
extracting correct FF values. 

The lower part of the diagram represents the doubly virtual Compton scattering (VVCS) off the proton. The 
amplitude of VVCS has two poles which are due to single proton in the intermediate state. The contribution of these 
poles to the box diagram is called the elastic contribution. Similar contributions appear from A(1232) and other 
resonances, Fig^ In the present paper we study the elastic contribution only (diagrams FigQJb) and (c), which we 
later on call box (in the narrow sense) and x-box diagrams). We also believe that the method derived here can be 
applied, with minor modifications, to the contributions of resonances. This will be the subject of a separate paper. 
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FIG. 1: Two photon exchange for the elastic ep scattering. 





In previously published papers the box diagram was calculated in several, more or less approximate ways. 

First is the soft photon approximation used in Tsai's paper (and also in Ref.[j|). It is assumed that the main 
contribution to the loop integral comes from the qi f=a 0, q 2 ~ Q and q\ ~ q, q 2 ~ regions. The reason is usually 
given that the integrand is strongly peaked at that points because of the infrared (IR) singularity. However, only 
the IR divergent part (which exactly cancels in the final answer), can be calculated in such a way. Instead, we are 
interested here in the IR finite part, and it is by no means obvious, that the main contribution to it comes from the 
same regions. 

The second way, which may be called "hard photon approximation" , assumes the most important region to be 
qi ?s q 2 fa q/2 This may be a reasonable assumption, however we believe it needs justification. In Ref.Q, after 
making the approximation the loop integral becomes divergent and an ultraviolet cut-off is inserted by hand. 

In the third group of papers the loop integral was calculated exactly, however using the special FF parameterization, 
monopole FF in Ref.[f| and the sum of monopoles in Ref.[||, 
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where rii and di are fitted constants. Thanks to this simple form, the resulting integrals are expressed through 4- and 
3-point functions, which were calculated numerically using computer program. 

However, such sort of calculation has the following problem. Since the FFs, entering the integral, are not known 
from the first principles, we should use some model for them or fit to experiment . The problem is that only 
t < region is accessible for the scattering experiments. While the t > 4M 2 region can, in principle, be studied in 
er +e + — ► p+p reaction, the unphysical region < t < 4M 2 is completely inaccessible. But the loop integral involves 
FFs in the time-like region (t > 0) as well. A natural question arises, how the (largely unknown) FF behaviour for 
t > influences the total result. In particular, parameterization used in Ref.0 does not even roughly reproduce FFs 
in the t > region. Indeed, unitarity requires that FFs have poles at t = m v where ray are masses of vector mesons, 
namely p-, to- and ip- mesons (neglecting their widths). But the numerical values given in Ref.0 are away from these 
masses. So, the fit is not suitable for t > and the results of Ref.|(| cannot be considered reliable before answering 
the above question. 

In the present paper we calculate the box and x-box diagrams in the most rigorous way. First, contrary to_Q,0, Q 
we evaluate all integrals exactly, without any restriction of the integration domain. Second, contrary to 0, Q we do 
not use any special assumptions on the FFs functional form. We perform the analytical integration to the maximal 
possible extent, resulting in 
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lCi J {t u t 2 )F 1 {t 1 )F J {t 2 )dt 1 dt 2 



(5) 



where K.ij(t\, t 2 ) is explicitly known. During the derivation of © we perform Wick rotation of the integration path, 
which is possible due to FF analyticity. As a result, the integration is done over negative t\ and t 2 only. Therefore 
we overcome the above-mentioned problem of finding out FF values at t > 0. Now the experimentally measured FFs 
can be used directly for the calculation of JSJ. Then we calculate two photon exchange amplitudes numerically using 
different FF paramcterizations and discuss the results. 



Of course, this is a vicious circle, since the accurate FF extraction from experiment implies taking into account the box diagram. 
However as a zeroth approximation we may neglect it; then some iterative procedure may be used to obtain precise result. 
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II. THE METHOD 



The box and x-box diagrams, as well as notation for particle momenta, are displayed in Fig|2| Thin line is for 
electron, thick is for proton. The proton and electron masses are M and m, respectively. We also denote 

P = (p + p')/2, K = (k + k')/2, q = p'-p. (6) 

The following relations arc fulfilled: 

qi = p" -p, q 2 = p" - p', k" =K±(P- p"), (7) 

upper sign is for the box, lower is for the x-box diagram. 
The amplitude corresponding to either of diagrams is 

_(a\ 2 f N(p")d 4 p" 

b?X ~ WJ J (q 2 - \ 2 ){q 2 - \ 2 )(k" 2 - m 2 ){p" 2 - M 2 ) ' { ' 



>t-bo 



where 

N{ P ") = u%(k" + m) lv u ■ U'T^q 2 )(p" + M)T v ( qi )U, (9) 

is defined by (0 and we use the notation a = a M 7 M . For the x-box diagram one should interchange 7„ with ~f v in 
Eq.®. The "photon mass" A is introduced in the denominator to avoid IR divergence. Though the electron mass is 
small compared to the characteristic energies involved, we will not neglect it in the denominator, first, for generality, 
second, since, as it is seen from the result, each of the diagrams diverges like lnm as m — > (but their sum does not). 

In general case the elastic ep scattering is described by six invariant amplitudes. However three of them are 
proportional to m, so in the m — > limit remain only three 0: 



M = -^T^u ■ U> (F l7 " - F 2 [^,^}^ + F 3 Kj^ ) U. (10) 



The invariant amplitudes Fi depend on two kinematic variables, v = APK and t = q 2 . In the first order of PT 
Fi(v,t) = Fi(t), F2{v,t) — F2(t), F3(v,t) = 0. The contribution of the box diagram to any of these amplitudes is 
given by the integral, similar to ©, but with another numerator, 



2 



N(p")= £ a^fmDfM), (11) 

i,j'=l 

where Aij(p") are some (rather cumbersome) explicitly known scalar polynomial functions ofp". To simplify the 
notation, from now on we drop the sum sign and the summation indices in the expressions like <|11[) and write them 
simply as A(p")F(q 2 )F(q 2 ). However the summation is always understood. 

There are at most four independent scalar functions ofp", say, p" 2 , pp" , p'p" and Kp" (the pseudoscalar combination 
^iiuarP'^P^p' K T cannot appear because of the parity considerations). The alternative, more convenient choice is 

p" 2 - M 2 , 

k" 2 - m 2 = tl+ * 2 ~* ± 2KP =f 2Kp", . 

h = q 2 = p" 2 + M 2 ~ 2pp" , ( 12 > 
t 2 = ql= p" 2 + M 2 - 2p'p", 




FIG. 2: Box diagram (left) and x-box diagram (right) with the proton in the intermediate state. 
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upper sign is for the box, lower is for the x-box diagram. Since Aij{p") are polynomials in p", they can also be written 
as polynomials in four variables (|12fl . 

Now it is easy to see that the integral (JHJ) can be reduced to the four basic integrals: 

AfrMmm^ (13) 



D 

where A is a polynomial, F(t) = F(t)/{t — A 2 ), and 
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= 1 
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= k" 2 ~ 
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= (fc" 2 - 


- m 2 )(p' 



(14) 

" 2 — M 2 ). 

If the maximal power of k" 2 — m 2 or p" 2 — M 2 in A{p") is greater than one, then other integrals can appear in the 
decomposition of (JSJ. However, they are expressed through l|13|) using symmetry considerations, see Appendix 1X1 

For the x-box diagram the integrals 7i, I 2 , I3, are the same as for the box diagram, but the integral I4 is not, since 
the relation between p" and k" is different. We denote the corresponding integral 1^. 

We will show here in detail the evaluation of the integral I4. Other integrals are evaluated similarly. At the end we 
write down the results for all of them. 

As it was discussed previously, the idea is to integrate over two of four integration variables, to obtain 

/„ = J }C n {t 1 ,t 2 )A(t 1 ,t 2 )F{t 1 )F{t 2 )dt 1 dt2 (15) 
with an explicitly known functions K, n . This expression can be used for numerical integration or further analysis. 

Evaluation of PC„ 

We will perform the calculation in the Breit frame. In this frame we have 



P = (iV4M^7, 0,0,0), ? = (0,0,0,,/=*), K = (g^fc, ^ 2 -%MlT~ t] > ^ ^ 
(the first is the time component, the following are x-, y-, and z- components). Let also 



p" = (£ + \^J^M 2 -t,pcos^,psm^,r]), (17) 
so d i p" = d£dr] pdp dcf>. In this notation 

t h 2 = e-p 2 -(v±v~t/2) 2 7 



p" 2 - M 2 = (£+±V4M 2 -t) -rf-p 2 -M 2 , (18) 



k" 2 -m 2 = ( -n - - - A 2 - p 2 - V 2 - m 2 - K 2 + 2pK x cos < 



We are to calculate the following integral: 

A(t 1 ,t 2 )F(ti)F(t 2 ) 



(p" 2 - M 2 ){k" 2 - TO 2 ) 
First we integrate over (j). The only quantity that does depend on <f> is fc" 2 . Using 



d£dr]pdpd(f). (19) 



1 f 2v dd> 1 



(20) 



2n J z- cos 4> ^/ z 2 - 1 
it is easy to verify that 
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FIG. 3: The integration path before (left) and after Wick rotation (right). 



In Ea. H20|) and further the following convention is used: we mean by v ' z 2 — a 2 the analytic function of z with the 
branch cut from —a to a such that V ' z 2 — a 2 > for z > a and \J z 2 — a 2 < for z < —a. Consequently, Ea. H21|) 
implies, that if is real, it has the same sign as the expression in the brackets under the radical. The integral 
becomes 

h = 2-K [ ^ Ut ti { ^} t2 \ didr l pdp= [ ^d^drjdp 2 . (22) 



(p" 2 -M 2 )D^i h p) 

Next we will integrate over £. Here we intend to use Wick rotation, so we need to study the analytic properties 
of the integrand. The FFs F(t) (and thus F(t)) are analytic functions of t everywhere in the physical sheet except 
the positive real axis. Namely, FFs have branch cuts running from t = (2m T ) 2 to +oo (m ff is the pion mass) 
@, and the additional pole at t = A 2 is due to the photon propagator; the physical value is F(t + iO). In the 
complex £ plane the corresponding singularities lie at £ 2 > p 2 + (|?/| — -y/^I/2) 2 . The factor l/(p" 2 — M 2 ) has 
two poles 1 and 2 at £ = — i\/4M 2 — t ± \J p 2 + if- + M 2 and has two branch cuts 3 and 4 running from 
£ = ^-j= T = ± y/(p— K x ) 2 + i] 2 +m 2 to ^-j=^= ± y/(p + K x ) 2 + i] 2 + m 2 . When integrating over £, we must pass 
by these singularities adding —iO to the masses in the usual way. The resulting integration path t is shown in FigEl 
on the left. 

Now we perform the Wick rotation to superpose the integration path with the imaginary axis. If the singularities 
1 and 3 lie at £ > and 2 and 4 lie at £ < 0, then the integration path can be rotated without crossing them, FigOU 
upper drawing. While singularities 2 and 3 always fulfill the above condition, it may not hold for 1 and 4- ln this 
case an extra terms appear when the path crosses the singularities, FigEl lower drawing. 

In general case wc may write 

/ = $ d£ + / A$ dt 9(00 ( i7 jfc - f) 0(-D 2 ) - (23) 

J£ J —too J 

-2m [ [{p" 2 - A/ 2 )$] 0(-{)0 (£ + |V4M 2 - t) S(p" 2 - M 2 ) d£, 



where A$ = $(£ — iO) — $(£ + iO) = 2<f>(£ — iO). It can be verified that during the integration according to l|23ll t\ 
and ti are always negative. For the first integral in the r.h.s. £ 2 < and thus t 12 — £ 2 — p 2 — [f] ± \/~-t/2) 2 < 0. 
As we will see later, the same is true for other integrals. This is a great advantage, since the FFs are well-known for 
t < (which corresponds to the elastic scattering), and much worse known for t > 0. 

Here we will change variables to make new independent variables (t l5 t 2 , £) instead of (p, r), £). It is important 
to take into account that p 2 > 0; for this purpose we insert 0(p 2 ) under the integral. This implies the condition: 

£ 2 > — | + tl 1 2 t2 — < ~ tl ^* 2 ' 1 = £ 2 , and the integration path along the imaginary axis becomes bounded. 
It is convenient to introduce the notation 

XM = \{^f-{\-^)t 1 t 2 , (24) 

™ _ 1 f t!+t 2 -t \ 2 (1 1 \. . 

Note that at t\,t 2 < one has x m > xm > %oo- 
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FIG. 4: New integration path £' . 



In new variables wc have 



p" 2 - M 2 = ^4M 2 - t + ti+^l = V4A/ 2 - t{£ - £ - C), (25) 



Dcj, = -V4m 2 -Wit -Ho) 2 -B 2 , (26) 

whcrc eo = tl+ ^- t (4m ._^ 4M2 _ t , S = feV^- C 1 = (l + 4^=7) • Therefore when $ is expressed 

through ti, t2 and £, it has, as a function of £, just one pole and one branch cut. It is worth noting that (|26[) may 
differ in sign from defined according to J2IJ- However along the integration path they are always equal. 
Rewriting 123fl in new variables and after long algebraic transformations (see Appendix [5]), we obtain 

dt ldt2 AfrMmm f e(Xoo) f g (27) 



2^~t V4m 2 - tV4A/ 2 - t \ y 'Je'(£-£o- C)^(£,-^) 2 -B 2 



*i,t 2 <0 



+2*fl(x m )e(-x 00 )fl(ti +t 2 -t) 



,-B (e-eo-C)^ 2 -^-^) 2 

1 



-2™ + 0(a;A/)e(-ioo)0(*i + *2 - *)] 



\JC 2 - B 



where i" is the integration path depicted in FigQJ 

The integration over £ now can be done analytically. The second integral, after substitution £ = £0 + B cos tp, is 
done with the help of Eq. H2U|) and gives 

* ^ - ^ (28) 



ho-B {£,-^-C)^B 2 -tt-t ) 2 Jo Scos^-C VC 2 - B 2 
The first integral is done using 

dz 1 , Cz-B 2 - VC 2 - B 2 Vz 2 - B 2 ... 

:ln 7^ „, , ^, „, = A ^)- ( 29 ) 



(z - C)Vz 2 - B 2 2VC 2 - B 2 Cz- B 2 + VC 2 ~B 2 Vz T ~B 2 

Here and below lnz is defined with the branch cut from z = to +00, so that for z > Imln(z + iO) = and 
Imln(z - iO) = 2tt. 

It is easy to show that the argument of logarithm in Ea. (|29|l can be real positive only if z is real, so all the 
discontinuities of A(z) lie on the real axis. Therefore 

A(£* - £0) - A(-£, - £0) - A(+oo + »0) + A(+oo - iO). 



Using 

v/(£* - Co) 2 - 5 2 = ^ (^fe - ^f=s) (30) 
after some simplifications we arrive at the final result of the form l|15|l with IC4 which is given below together with 
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other K, 7l 



JC 2 {tiM) 



0(Xoo), 



2^-t{Am 2 -t) 1 V ' V 



27ri0(a; m )0(-a; oo )0(ti + t 2 - t) 



2 v /-t(4Af 2 -t) 



(x„c) In 



fl(a?oo) in (e + 2 ^"-t ) |1* ~ 2 ™ QMoi-x^eih +t 2 -t)\, 



2JR+ 



27ri ^(xoo) + 2™ Oi-XooMh +t 2 - t)[6(x m ) + 8{x M )]} , 



2,/RZ 



(Xqo) In 



+ 2™ + t 2 - t)[6»(x m )6»(-a;M) + 20{R-)0{x M )0(y + t - AM 2 )}} , 
where = iy/Xoc and 

fl ± = ( «i+^-« ) 2 ((„ T t) 2 - 16m 2 A/ 2 ) - ^(i/ 2 - (4m 2 - <)(4Af 2 - <)), 
V /r R± = o'l-Ril 1 ^ 2 ; where a = -sign(ii + t 2 - t) if i?± > and a = -i if R± < 0. 



(31) 
(32) 

(33) 

(34) 
(35) 
(36) 



III. NUMERICAL 



It is convenient to distinguish three integration regions: i) Xqc > ii) xm > > Xoo and iii) x m > > xm- In 
general case the integrals I 2 , I4 and Ii x contain logarithmic terms like In ^ originating from the integration over 
x m > > xm region. However in actual calculations after adding box and x-box diagrams these logarithms always 
cancel; it can be shown to be the consequence of gauge invariance. Therefore we may put m = if the amplitude is 
evaluated as a whole (box + x-box). 

Before calculating ep — > ep amplitudes, the following cross-check can be performed. If we omit A(ti,t 2 ) and set 
F(t) = 1 in Ea. l|13|) then the integrals /„ can also be done by usual Feynman parameters method. We have calculated 
them numerically with different values of M , m, A, t, v, and made sure that both methods give identical results. 

Now we turn to evaluation of invariant amplitudes F 2 , F3. Each of them looks like 



Fi = ai]n\ + bi + o{\). 



(37) 



The first (IR divergent) term does not contribute to the cross-section if the radiation of soft photons is taken into 
account. To extract both IR divergent and IR finite part and to simplify the calculation we used the following 
procedure. We calculated the integrals needed at several small but non-zero A and fitted obtained values with 3- 
parameter fit: 



Fi = cii In A + hi + CiX. 

Though the third term vanishes as A — ► 0, it turns out to be necessary to obtain accurate results. 
The results obtained by Tsai 01 are 



TT 



0, 



(38) 



(39) 



where K(pi,pj) = (jpi ■ Pj) J - 



p y = piy +pj(l — y); the superscript (T) hereafter indicates Tsai's result. In 



spite of the approximate nature of this result, the terms proportional to In A are exact; they are 



a l,2 — a l 2 



2a , v - t 
— F 1 2 In . 

7T ' V + t 



(40) 



Our numerical calculation has given the same results for these terms. 
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FIG. 5: Two photon exchange corrections SGm/Gm, 5Ge/Gm and F27 as indicated on the figures, for Q 2 (a) 1 GeV 2 , (b) 3 
GeV 2 , (c) 6 GeV 2 , (d) 10 GeV 2 , using form factor parameterization: from Ref.0 (solid), from Ref. llOl (dashed) and dipole fit 
(dotted). 



Instead of F\ t % wc introduce linear combinations Gm and Ge exactly as this is done for FFs, Ge = F\ + 4-^2-^2 
and Gm = A + F 2 . To see the relative size of corrections with respect to the Born approximation, we consider the 
following quantities: 

SGm/Gm = (Gm — G^p — Gm)/Gm, 

SGe/Gm = (G E -G ( P -G e )/Gm, (41) 
Y 2l = (v/4M 2 )(F 3 /G M ). 

The results of Tsai are subtracted, since they were used (the therefore already taken into account) in the analysis 
of almost all experimental data. This is not the same as just drop the logarithmic term, since l|39|l contains non- 
logarithmic terms as well. 

Note that in contrast to Ref.|(|, we normalize SGe by Gm, rather than Ge- This is done because of great discrepancy 
in Ge as extracted by two methods; in addition, as extracted by polarization transfer method is close to zero at 
Q 2 f=a 6 GeV 2 and SGe/Ge becomes large though SGe itself is small. On the other hand, the uncertainty in Gm is 
much smaller. 

We have calculated the two photon exchange amplitudes at Q 2 = 1,3,6 and 10 GeV 2 using three FF parameteriza- 
tions: 

• dipole fit 

• fit from Ref.0 (under assumption Ge/Gm = const) 

• fit from Ref.[Tol| (to data obtained by polarization transfer) 

The results are plotted in Fig|3 versus customary parameter e (virtual photon polarization), e = . 

We see that in general results agree with those obtained in Ref.[(|. This is not surprising, because in Sec ITU wc 
actually prove that given the FF parameterization which is good for space-like region (t < 0) and has correct analytic 
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properties, the result will not depend on the quality of the fit for t > (a priori this was not clear). In some sense we 
justify the approach of Ref.Q. Nevertheless, our approach is more general, since we can use any FF parameterization. 

We also see that the dependence on the choice of FFs is small for Q 2 < 6 GeV 2 , however for the generalized electric 
FF at Q 2 = 10 GeV 2 the difference reaches 50% (for small e). 

The absolute value of corrections increases with Q 2 up to ~ 3% for SGm/Gm at Q 2 = 10 GeV 2 . 

The main purpose of this paper is to present new approach for evaluation of two photon exchange diagram. That is 
why we will not analyze the effect of these corrections on the FF measurements. Such analysis requires evaluation of 
bremsstrahlung corrections as well, and thus detailed consideration of experimental conditions is needed. We postpone 
it to another paper. 

IV. SUMMARY 

We study the two photon exchange in the clastic ep scattering. In this paper we calculate box and crossed-box 
diagrams with proton in the intermediate state (elastic contribution). Our approach has two main advantages: 

• the amplitude is expressed via proton form factors in the space-like region only, Eq.(|SJ). The previous calculations 
required form factors in the time-like region as well, where they are poorly known. 

• any suitable form factor parameterization in the space-like region can be used to evaluate the loop integral. 
In previous calculations only particular kind of parameterization could be used, since the loop integral was 
expressed via 4-point functions. 

Similar approach can be used to evaluate A(1232) and other resonances contributions, which is currently under 
investigation. At present only the A contribution was considered in the literature using approach of Refs.0,0 
and simple dipole parameterization of the N — * A transition form factors. 

We have calculated two photon exchange amplitudes using form factors as extracted by Rosenbluth and by polar- 
ization transfer methods. The dependence of the result on the choice of form factors is small except the generalized 
electric form factor for Q 2 > 6 GeV 2 . 

The same method can be applied to the scattering off any spin-1/2 particle, such as neutron or 3 He. 

APPENDIX A 

In this Appendix we show how some integrals, that can appear in the decomposition of JHJ, arc reduced to four 
basic types (fTTfl) . 

Consider the integral 

I A( tl ,t 2 )^-dY = f A(t u t 2 ) ( p+ l7, k Z- M2 <PP" = (Al) 
= J Aih^yi^^p" + J A{ti,t 2 )d i p" - 

-2(P + KY J Ait^yS^d^p" . 

The first and the second integral in the r.h.s. already have required form. Since t± >2 = (K ± q/2 — k") 2 , the last 
integral depends only on g M and K^, thus 

/ A(ti,t 2 ) W 7^^d i p" = o% + PK^, (A2) 

where 

a=±jA(t u t 2 ) 1 ^d* P >> P=^ t f A(h,t 2 ) [ A ^-J 2 + 2 ) d A p", (A3) 

so finally 

fA( tl ,t 2 )^EBd^" = | (l+ 5;^) fA( tl ,t 2 )j^d^" - ^fAit^W- (A4) 
Similarly one obtains 

f A(h,t 2 )^E^d*p" = i (l + jj^j) j A(t u t 2 )p^d^ P " -^-jAQuhW (A5) 

and 

/ A{t x ,t 2 ){k" 2 - m 2 )d 4 p" = J A(t u t 2 )(p" 2 - M 2 )d 4 p" = jA(t 1 ,t 2 ) t -i±^d 4 p". (A6) 



FIG. 6: Regions in (ti,t2) plane, where x m > (hatched). 



APPENDIX B 



In this Appendix we show in detail the derivation of Eq. (j27(l . 

Let us transform two last integrals in Ea. l23fl . It follows from Eq.J2U that if D 2 < then always p 2 > 0, so 
9(p 2 )9(-Dl) = 9(-D 2 ), and also 

Dll = (P 2 + m 2 + rf - K 2 x f + 4K 2 (m 2 + if) > 0, (Bl) 



2y/4M 2 -t 

so 9 ( 2 ^ 4 ^ f 2 i ~ £j is constant (cither 1 or 0) when < 0. Therefore 



^ 9 ^ (iVm^t - ^) 0(B 2 - - io) 2 ) = (B2) 
= 9(£)9(4m 2 -t x - t 2 )9(x m )9(B 2 - (£ - £o) 2 )- 



(x m ) can be inserted since < B 2 — (£ — £ ) 2 implies < B 2 = ( ^i^ t ) x m . Rewriting x m as 



4m 2 -r m ~~ 4m 



bl(^ 1± ^) 2 + H^) 2 -^ 2 (B3) 



% ' -x. 



it is easy to see that x m is positive in two regions, bounded by hyperbola x m = (FigEJ) . The condition 4m 2 — 1\ —t 2 > 
selects lower- left region, where ti,t 2 < 0. Thus 

e(Am 2 -tx- t 2 )8{ Xm ) = e(-h)e(-t 2 )e(x m ) (B4) 

and 

°(P 2 )°(00 (j^mrn - e) 0{-D 2 ) = 9(-tx)9(-t 2 )9(x m )9(09(B 2 - (f - ^o) 2 ). (B5) 
To transform the last integral, we write 

O(p 2 )0(-£)O {£ + y^M 2 - t) 5{p" 2 - M 2 ) = (B6) 
= 0(e + Xoc )0(-OO (e + yui 2 - t) 5 (i^fAM 2 - t + *ii^=*) = 

= e(x M )0(tx +t 2 - t)9(UI 2 -tx- h) 7 ^ w =6 (t + Jfife) = 

= tfHiw-w^wti +t 2 - t) 7 ^ T = t s (e + ^f^) . 

The last equality was obtained using Ea. (|B4(l with m replaced by M. Again we have tx,t% < 0. 
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After that the integral I4 becomes 

f + ZOO 



f dt dt ( r+ % °° f 



*1,*2<0 

2?ri 



0(ti +h- t)6{x M ) f [{p" 2 - A/ 2 )$] 8 + #} . 



Also 



and finally 



(B8) 



V4:M 2 - t 

Now we introduce the integration path £' that passes by the singularities on the right, Fig0| After that 
f 0(? + xoo) <5>d£ = 6( Xoo ) I jT $ d£ - J A$ 0(O^(S 2 - (£ - Co) 2 )- 

-2™ 0(£ o + C) J [(£ - £ - C)$] - - C) d£ 
Combining this with the previous equation, we have 

r fl(xoo) / *de + fl(a: m )fl(-Xoo) / A$0(O0(B 2 -(£-£o)*K- 



ti,i 2 <0 

-2 7 r l [0(x oo )0(i - *! - 1 2 ) + e{x M )e{t 1 + 1 2 - 1)\ J [(e - & - c)$] <y(c - & - co d$ 

To proceed further, we note that 

B 2 (C - £o) 2 | £=0 = (4m ,_ t ^-t) ^ - 3^x m < (BIO) 
for Xao < and x m > 0. Therefore 

e(x m )6(- Xoo )9(06(B 2 - (£- e ) 2 ) = (Bll) 

= Oix^ei-x^eiZoWB 2 - (f - Co) 2 ) = 

= Oix^Oi-x^Oih + f a - t)^(B 2 - (£ - Co) 2 )- 

6{ Xoo )6{t -h- t 2 ) + e{x M )e{t 1 +t 2 -t)= (B12) 

= 0(Xoo) + (0(.T M ) - 0(Xoo))»(*l + *2 - *) = 

= ^Gtoo) + e{x M )e{-x 00 )e{t 1 + t 2 -t) 

h= J j^xoo) jf * # + ofrmm-xooWh + * a - t) J a$ ^ - (e - £o) 2 )d£- (Bi3) 

ti,t 2 <0 

-2 7 ri[^( : r 00 ) + 0(ato)0(-a:oo)0(*i + 1 2 - *)] / [(£ - f - C)$] <5(£ - £ - C) d£ 
Substituting, according to (|2*2*|l 

V^^WUP - 1 (£-£o-CV(£-£o) 2 -£? 2 ' 
y/totty/m* - 1 (£-£ -CVfl 2 -(£-&) 2 ' 

into (iBTBTl . we obtain Ea. (|2Tjl . 
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